Slow relaxation of conductance of amorphous hopping insulators 
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We discuss memory effects in the conductance of hopping insulators due to slow rearrangements of 
structural defects leading to formation of polarons close to the electron hopping states. An abrupt 
change in the gate voltage and corresponding shift of the chemical potential change populations 
of the hopping sites, which then slowly relax due to rearrangements of structural defects reducing 
the density of states. As a result, the density of hopping states becomes time dependent on a scale 
relevant to rearrangement of the structural defects leading to the excess time dependent conductivity. 

PACS numbers: 73.23.-b 72.70.+m 71.55.Jv 73.61. Jc 73.50.-h 73.50.Td 



I. INTRODUCTION 

A peculiar transport memory effect has been observed 
in many hopping insulator^i^i^i^i^: Shifting the system 
either from the equilibrium or from its steady transport 
state produced by, e.g., a sudden change of the gate 
voltage, 5Vg^ increases the conductance a of the sys- 
tem, the effect not depending on the sign of the change. 
The cF{5Vg) dependence shows a characteristic memory 
cusp (see Ref. l6| for experimental details and a review) 
which may persist for a long time. Explanations of this 
appealing observation can be grouped into two major 
classes which are often referred to as intrinsic and ex- 
trinsic mechanisms. The former one attributes mem- 
ory effects to slow dynamics of strongly correlated elec- 
trons subject to quenched disorder and thus forming a 
Coulomb glass^^^. The connection between the features 
of a glassy phase and the formation of a Coulomb gap 
was revealed in Ref. HO, suggesting that within the loca- 
tor approximation the correlated electronic system maps 
to the Sherrington-Kirkpatrick spin glass. A direct re- 
lation between the long range electron correlations and 
the formation of the exponential distribution of deep 
electron states which is characteristic to glassy systems 
was also demonstrated in Ref. [ll|- In a Coulomb glass, 
memory effects reflect delayed formation of the Coulomb 
gap in the Efros-Shklovskii (ES) regime of the variable 
range hopping ( VRH) Recently we have demonstrated 
how "slow" many-electron fluctuators can be formed in a 
Coulomb system and analyzed their influence on the 1// 
noiseJ^ Such fluctuators can also lead to slow relaxation 
and memory effects. 

An extrinsic scenario assumes that electronic memory 
effects are caused by slowly relaxing atomic configura- 
tions influencing conducting channels and was first pro- 
posed in Ref. n4 to explain the G{Vg)-cusp in granular 
Au films. A possibility that polaron effects may be re- 
sponsible for slow relaxation in hopping conductors was 



also discussed in Refs. |5|Ui 

The Coulomb glass is formed well below the tempera- 
ture To entering the ES law for the VRH as 



a = (Joe 



-(To/T)i/2 



(1) 



Experiments^ has shown that in this regime the 1// noise 
intensity is strongly correlated with To proving electronic 
nature of the low- frequency noise. However, the conduc- 
tance memory cusps and their relaxation (ageing) where 
observed in the systems which do not exhibit ES behav- 
ior, but are apparently subject to structural disorder— lii. 
Various memory effects were observed in metallic gran- 
ular structures-i^ii^, which also possess a high degree 
of structural disorder. This calls for careful examining, 
whether the memory effects and the "two-cusps" G{Vg)- 
dependence^ can be explained as a result of slow relax- 
ation in the structurally disordered atomic matrix. 

In this Letter we present a model that may serve as a 
step toward a quantitative description of memory cusps 
based on an extrinsic mechanism. We show that due 
to slow relaxation of atomic structure, polaron clouds, 
which form near the hopping sites, suppress the bare 
electron density of states (DOS). Changing gate voltage 
shifts chemical potential, removes the polaron screening, 
and, thus, increases the hopping conductivity. As the 
atomic structure adjusts itself with time to the new po- 
sition of chemical potential, the conductivity relaxes to 
its quasi-stationary magnitude. 

Atomic structural relaxation is described by two level 
systems (TLS)^^ The TLS model successfully de- 
scribes thermodynamics and kinetics of amorphous solids 
at low temperature. It suggests that there exist atoms 
or groups of atoms undergoing tunneling motion and 
characterized by the broad universal distribution of their 
parameters. Since all materials where the memory ef- 
fects in conductivity were observed are strongly disor- 
dered one would expect that two level systems (TLS) 
should exist there similarly to other glasses and disor- 
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dered materials These two-level systems interact with 
conducting electrons because they possess the dipole 
moment. In this manuscript we examine the effect of 
electron-TLS interaction on the non-equilibrium conduc- 
tivity. We show that this interaction results in the non- 
equilibrium behavior of conductivity which is qualita- 
tively equivalent to the experimental observations, i. e. 
increase in conductivity after gate voltage application 
with its subsequent logarithmic relaxation to the equi- 
librium value. This theory uses the previous work^i^ 
where the similar non-equilibrium behavior of the dielec- 
tric constant in amorphous solids^^'^ has been explained 
using the TLS interaction. Since TLS parameters are 
quite universal from material to materiaU^ we can use 
these parameters for quantitative estimates which show 
that our theoretical predictions are consistent with exist- 
ing experimental data. 

The paper is organized as following. In Sec. [Ill we ac- 
curately define the conditions where our consideration is 
applicable, i. e. electron-electron interaction can be ne- 
glected, while electron-TLS interaction is significant. In 
Sec. mil the non-equilibrium behavior of electron density 
of states caused by their interaction with TLS is derived. 
In Sec. IIVI the non-equilibrium behavior of conductivity 
is described and compared with the experimental data. 
The results of the manuscript are summarized in conclu- 
sive Sec. |Vl 



II. WHEN THE ELECTRON-ELECTRON 
INTERACTION CAN BE NEGLECTED? 

In this manuscript we ignore the effect of electron- 
electron interaction on the non-equilibrium behavior of 
conductivity. This is possible only under specific condi- 
tions when the electronic interaction is weak compared 
to their characteristic energies and there is no slow relax- 
ation within the electronic subsystem. According to var- 
ious considerationsi2J^ the slow relaxation in electronic 
subsystem can take place only below some critical tem- 
perature Tg which is defined as the electronic glass tran- 
sition temperature in Ref. [T2|. This transition tempera- 
ture is defined by the Coulomb gap energy 

UbTg ^ Ac. (2) 

Below we set ks = 1. 

The Coulomb gap depends on the electronic density of 
states go^ dielectric constant system dimension dand 
electron localization radius a. We believe that the dimen- 
sionless parameter x = g^a^'^e^ /n = g^a^T^ is small, i. 
e. 

X = ^oa''To<l, (3) 

where the characteristic temperature Tq in the Efros- 
Shlovskii hopping conductivity Eq. ([T]) is given by 

kBT^ . (4) 



Then the Coulomb gap is defined as 

Ac=Tox^. (5) 

If temperature exceeds the Coulomb gap energy then 
electronic relaxation is fast. It is also useful to notice that 
for the electronic excitations with typical energy E their 
characteristic interaction Ue ^ e^(^o^)^ is smaller than 
the energy E until E > Ac and the condition Ue ^ E 
serves as the definition of the Coulomb gap. So when 
T > Ac the interaction of representative excitations with 
energy E >T is less than their energy so it can be treated 
as weak. 

If the system is near metal dielectric transition^ and 
the electron localization radius is large, then one can 
possibly have the opposite limit x — g^CL^~^e^ > 1, 
This is the case which takes place in silicon Mosfets in- 
vestigated in Ref. [IB. Under those conditions the only 
Efros-Shklovskii variable range hopping law Eq. ([T|) is 
observed. We believe this is not the case for the sys- 
tems of interest ^i^i^i^i^ where the different conductivity 
behavior is observed more close to the Mott's variable- 
range hopping law, which suggests x < 1-^^ 

Therefore it is not possible to interpret the memory ef- 
fects in conductivity at KbT > Ac without the involve- 
ment of electronic interaction with the extrinsic slowly 
relaxing defects like TLS. Our theoretical study is re- 
stricted to this "high temperature" situation. Since the 
temperature Tq can change within the range of 1 — lOOK— 
our assumption that the experimental temperature 4K 
exceeds the Coulomb gap energy Eq. ([5]) does not con- 
flict with the common sense. We believe that at least 
some of experiment^i^i^i^i^ are performed in this temper- 
ature range. Unfortunately, existing experimental data 
do not permit us to answer the question whether experi- 
mental temperature is above or below the Coulomb gap. 
Although the absence of the Efros-Shklovskii conductiv- 
ity temperature dependence Eq. agrees with our as- 
sumption that the electronic interaction is not important 
this does not prove the absence of Coulomb gap and elec- 
tronic glassy state.^^. We therefore suggest additional 
experiments which can help to investigate the Coulomb 
gap in the system which is the necessary prerequisite of 
the Coulomb glass state. Such experiment can be made 
for instance using scanning tunneling microscopy. If the 
Coulomb gap in the density of electronic states will not be 
observed at experimental temperatures then one should 
expect that the extrinsic mechanism is in charge for the 
memory behavior in conductivity. 

HI. ELECTRON DENSITY OF STATES 
AFFECTED BY ELECTRON-TLS INTERACTION 

Our consideration is based on the concept of 'two level 
systems' (TLS)^ taking as its central hypothesis the as- 
sumption that in a system with quenched disorder a cer- 
tain number of atoms (or groups of atoms) can occupy 
one of (at least) two equilibrium positions. These atoms 
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therefore move in a double- well potential created by their 
environment and characterized by the asymmetry en- 
ergy (difference in energy minima) and by the height and 
the width of separating barrier. The atoms comprising 
the TLS change their configuration either by tunneling 
through- or by thermally activated hopping over the bar- 
rier. Randomness in the heights and widths of the TLS 
barriers gives rise to exponentially broad distribution of 
structural relaxation times. 

A localized electric charge polarizes the neighboring 
TLSs, which thus acquire the electric dipole moment, /l^, 
and screen the original bare charge, see Fig. [1^. The 
interaction energy of TLS with the localized charge is 
estimated as 

U{t) ^ e{iJ.-T)/K.r^ , r:^a, (6) 

where r is the vector connecting the center of this TLS 
with the position of localized electron state, a is the local- 
ization length, and is the material dielectric constant. 
Polarized TLSs form a polaron cloud around the localized 
electron and create a polaron gap for electronic excita- 
tions, implying that the states with the single-electron 
energy (j) < U{ri) cannot be excited any more. 

In what follows we will find the correction to the elec- 
tronic DOS, go^ due to electron-TLS interactions in the 
lowest approximation in the TLSs density. To this end 
we first determine the change in DOS due to a single 
neighboring TLS characterized by the given energy split- 
ting, and the relaxation time, r. The target correc- 
tion to the electron DOS is then obtained by averaging of 
the single TLS contribution with respect to all possible 
neighboring TLS positions and E and r. 

Consider the correction to the density of electronic 



respectively. The time-dependent correction to the den- 
sity of states taken at the certain time t can be associated 
only with those two level systems which have the relax- 
ation time r longer than the time of the experiment t. 
That time must also have upper restriction by some maxi- 
mum TLS relaxation time Tmax which serves as the upper 
cutoff in the TLS logarithmically uniform distribution 
over their relaxation times P(A,r) = PoQ{rmax — 't)/t. 
Within the logarithmic accuracy one can account for the 
above constraints introducing integrated time dependent 



states with some energy e caused by electron-TLS in- 
teraction. Since the concentration of TLS is small we 
can assume that only one most closely located two-level 
system is important, while the probability that two TLS 
are significant is much smaller than unity. Then we con- 
sider electron interacting with the single TLS using the 
standard Hamiltonian 

Hpair = (t)n^ UnCF^ + Acr^, (7) 

where n = 0, 1 is the electron population operator, the 
spin 1/2 z-projection operator cj^ = ±1/2 describes two 
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FIG. 1: a - Polaron formed by an electron site and an adjacent 
TLS; b - Transitions of the TLS between its states shifts the 
electron energy with respect to the Fermi level and can bring 
it outside the hopping band having the width Sh- 

states of the TLS, (j) is the electronic energy in some lo- 
calized state, A is the energy of the two level system and 
U is the charge-dipole interaction of the electron with 
the TLS. In the thermal equilibrium the electron excita- 
tion energy e = (j) -\- Ua can take values (j) ±U/2 with 
the probabilities defined by the equilibrium Boltzmann 
factors 



(8) 



I 

TLS density as 

P{t) = PoHt/Tm). (9) 

This expression clearly demonstrates the nature of loga- 
rithmic time dependence of TLS contribution to the elec- 
tronic of states similarly to previous work^^. 

The correction to the density of states of the particu- 
lar electron with energy e caused by its interaction with 
neighboring TLSs can be expressed as 



exp (-^) + exp ^ ^/'^ u/2 ^ 

exp (^) + exp (z A/2-^/2-0 -| ^ ^ A^^ ^ AZ2+^^' 

exp (^) + exp ^^/2±^y2^^ 

exp (^) + exp ( rN2-u/2-, ^ ^ ^ ^N2±m^^^ ' 
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(10) 



r 



where the sum is taken over ah pairs made of an elec- 
tron i and a TLS j characterized by energies and Aj, 
respectively, Uij is their interaction and V is the system 
volume. Probabilities are defined using Eq. (|8|) with 



J 



the substitution = 0^, A = Aj and U = Uij. The 
summation in Eq. (p!Q|) can be replaced with the integra- 
tion over electron and TLS energies and the distance R 
between electron and TLS. This yields 



d'^R / # / dA 

J — OO J — oo 



(1 + e-i')5{€ -(j)- U/2) (1 + e-^)5{e - + U /2) 



■ e T 



A e-U . 

■ eT (1^ e T j 



5{e-4>) 



(11) 



where qq is the electronic density of states at energy e 
taken in the absence of interaction with TLS and d is the 
system dimension. Remember that U = e(/i, Yi)/{nR^) is 
the dipole charge interaction between electron and TLS 
possessing the dipole moment /i and averaging is per- 
formed with respect to random directions of a TLS dipole 
moment. After integration over electronic energy (j) and 
TLS energy A in Eq. (pTj) we obtained 



j d^R ^In 



< Sg{e) >= goP{t)T 

cosh(g/T) + 1 y 



cosh{e/T) + cosl: 



(12) 



The most important electrons are those contributing to 
the hopping conductivity. These electrons have energy 
e of order of the hopping energy Sh which always ex- 
ceeds the thermal energy. Therefore assuming e > T 
one can approximate the logarithm under the integral 
in leading order in 1/T as In cosh(£/T)+i 



cosh(£/T)+cosh(C//T) 



11^1-^1/ 



rp -0{\U\ — e). The final expression for the correc- 
tion is similar to the one used in earlier work o^^i^^ for 
the Coulomb gap and the dipole gap in the density of 
TLS states. However our derivation of the correction is 
more general because our expression Eq. (p!2|) can be used 
at arbitrarily temperature while the earlier derivation is 
valid only in the low temperature limit. Finally the cor- 
rection to electronic density of states can be expressed 
as 



-2goPoC J dr [U{r) - e] 0[U{r) - e] 



5g{e,t) 



where /i is the typical dipole moment of a TLS and Tm 
is some characteristic minimum time associated with the 



gate voltage application. Only the electrons with the 
energy e ^ Sh = T'^l^a~'^l^g~^l^ influence the hopping 
conductance. We believe that at £ ~ e/^ ~ 30K Eq. (p!3|) 
is still applicable. 



IV. ANALYSIS OF EXPERIMENTAL DATA 

Equation ([13]) can be used to interpret the experimen- 
tal data only if all relevant TLSs were initially out of 
equilibrium. In other words, all electrons contributing to 
hopping must be surrounded by non-equilibrium TLSs. 
The equilibrium can be broken due to the application of 
the gate voltage. This can be realized when all electrons 
contributing to hopping are "newcomers", i.e., they are 
brought to the vicinity of the Fermi level by application 
of the gate voltage, Vg. For that the shift of the Fermi en- 
ergy by the gate voltage must exceed the hopping energy 
Eh' Then those electrons break the equilibrium in their 
neighboring TLS randomly changing their energies by the 
scale of their interaction with those TLS. All relevant 
TLS with the energy of order of e/^ coupled to entering 
electrons (or holes) by the interaction having the same 
order of magnitude experience the jump in their energy 
induced by entering electrons. This moves them all out of 
equilibrium. Their relaxation leads to the polaron shift 
of electron energy out of the Fermi energy thus reducing 
the conductivity. Thus the condition IS^Ep ^ Sh defines 
the width of the cusp in the non-equilibrium conductiv- 
ity as a function of the gate voltage. One can show using 
Ref.0 that at those conditions the shift of TLS energy in- 
duced "directly" by the gate voltage is still smaller than 
because of the small TLS dipole moment /i ~ 2 • 10 
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g^gi/2^^3/219 compare this estimate with the 

experiment extracting eh from the expression for the con- 
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ductivity, a ~ ctq exp(— s/^/T). Let us put ctq equal to the 
minimal metallic conductivity (~ 10~^ Ohm"^) and use 
available experimental results of Ref. H for the represen- 
tative sample with resistance of 3.8 MOhm at T = 4.1 K. 
The assumption about ctq is justified by the experimental 
observations (see Ref.^^ and references therein). 

Then the hopping energy is Sh = Tln(cro/cr) ^ 2.1 - 
10~^ eV. This value agrees qualitatively with the shift of 
the Fermi energy 3 — 5 meV associated with the value of 
the gate voltage Vg at which the memory cusp is affected 
as was estimated earlier in Ref. 0. Thus the suggested 
mechanism agrees with the experimental observations. 

The correction to the conductivity can be estimated as 
[{e/T)Sg{e,t)/go]^^s^. In this way we get 



6a{SVg,t) Poep 



1/2 



In I — 



In 



^0 

a{Ty 



(14) 



The time-dependent factor C = ln(t/rm) contains the 
measurement time as t and the inverse sweep rate, Tm ^ 
Th = Toe*^^^/^^^^^, and affects only the amplitude of the 
peak, but not its shape. Here is the characteristic time 
of the variable range hopping and the preexponential fac- 
tor To is of order of Ips at experimental temperatures. 
This can explain the experimentally observed indepen- 
dence of the dip shape on the sweep rate^. 

Similar considerations apply to the ES VRH with the 
proper renormalization of the hopping parameters. Now 
the typical energy scale optimizing the hopping rate is 
Sh = (ror)i/2, where To = /Se^Ka, th = Toe(^o/^)''', 
and ~ 2.8 is a numerical factor—. 

Let us now discuss the available experimental data in 
light of the above theory. According to our previous es- 
timates, £h ^ 2.1 ■ 10~^ eV. The typical dipole moment 
can be estimated as /i ~ 2 • 10~^^ erg^/^cm^/^, which is 
close to typical dipole moments of TLSs in glasse a-*^^!^^ . 

We are not aware of independent measurements of the 
TLS density, Pq, in the materials under consideration. 
In principle, this quantity can be determined, e.g., by 
measurement of the low frequency dielectric constant at 
T < 1 K, where it should depend on the temperature 
logarithmically, the slope of the logarithmic dependence 
being P^ji^ / n. This dimensionless quantity turns out 
to be almost the same, ~ 0.3 • 10~^, in many materi- 
als with strong quenched disorder, ^e^^^^ and ref- 
erences therein. Apparent universality of this quantity 
in such materials was attributed to the interaction be- 
tween TLSs^iiS^. Assuming that the material studied in^ 
belongs to the same "universality class" as the materials 
with strong quenched disorder we estimate Pq as 0.6-10^^ 
erg~^cm~^, as in the most oxide glasses^^. We set ti ^ 1^ 
following Ref. 9. 

Now we can check our assumptions regarding 3D ar- 
rangement of TLSs forming the polarons and regarding 
lowest approximation in the TLS density. Estimating the 
polaron radius, f, as as the length at which the electron- 
TLS interaction, e/i/^^f^ is comparable with the typical 
electron energy we get f ^ 1 nm that is much less 



than the sample thickness I. The average number of the 
TLS forming polarons, i.e., located within the polaron 
radius and having E < Sh^ M ^ (47r/3)Po£:/if^ turns out 
to be ~ 10~^, i.e., much less than 1. Thus the lowest 
approximation in the TLS density is valid. 

Using the above estimates and Eq. (p!4|) we get 
Sa{t)/a ^ 0.02 ln(t/rm). Thus we predict the logarithmic 
relaxation rate of conductivity r = dln{a)/dln{t) ^ 0.02. 
According to the experimental daid^^ at T = 4. IK and 
for the sample thickness equal to lOnm the conductivity 
changes by about 8% during two decades in time so we 
can estimate it as r ^ 0.015. Thus our theory agrees 
with experimental data reasonably well. 



V. DISCUSSION AND CONCLUSION. 

We have presented a simple model of slow dynamics 
of hopping conductance in structurally disordered hop- 
ping insulators. It takes into account rearrangements 
of the dynamic structural defects, TLSs, leading to for- 
mation of polarons close to the electron hopping states. 
The model qualitatively explains both the logarithmic 
relaxation and memory effects, and provides quantita- 
tively reasonable estimates of the time-dependent non- 
equilibrium change in conductivity, S(j{Vg^t) (see, e.g.^. 
Fig. 2). 

The dependencies of this quantity on different parame- 
ters - electron concentration, controlled by the gate volt- 
age; magnetic field, ^, and various protocols of breaking 
down the system equilibrium - are encoded in the log- 
arithmic factors In (cr/cTo) and ln(t/rm) while the tem- 
perature dependence enters as a power law through the 
energy Sh- Sa (x ' . That leads to the main tempera- 
ture dependence oc T~'^l^ and oc T~^/^ for the Mott and 
ES VRH, respectively. Thus our theory explains the fast 
increase of the non-equilibrium raise of conductivity with 
decreasing the temperature. 

The following note is in order. We demonstrated 
that "slow" excitations induced by structural disorder 
can indeed be responsible for the double-dip memory 
effect in hopping semiconductors in analogy to those 
in glasses^iiS^. The effect we considered is due to po- 
larons formed by the structural excitations. Yet, one 
should bear in mind that a similar polaron effect ap- 
pears in pure electronic models. In particular, the po- 
larons formed from pair excitations were considered in^i. 
In its turn, "electronic polarons" can be formed from 
slow relaxing electronic "aggregates" discussed in Ref. 
il3i. The presence of structural TLS in InO films is indi- 
rectly supported by the fact that the low-frequency noise 
in these materials does not increase significantly under 
the disorder-driven metal-to-insulator transition. 

Further experimental verification is necessary to de- 
cide whether the non-equilibrium behavior is associated 
with the extrinsic or intrinsic model. A regular approach 
to attain this goal should be based on techniques affect- 
ing differently the structural and electronic degrees of 
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freedom. One of the possibilities is measurement of the 
AC Unear response such as simultaneous measurements 
of the attenuation and velocity of the acoustic waves in 
transverse magnetic field. The expected effect is caused 
by transitions in structural two- well configurations, while 
the electronic transitions are strongly suppressed by the 
magnetic field^^. Alternative experimental verification 
can be performed using canning tunneling microscopy^ 
as described in Sec. HIl 
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